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ABSTRACT 

We report on third epoch VLBI observations of the radio-bright supernova SN 201 ldh 
located in the nearby (7.8 Mpc) galaxy M51. The observations took place at t = 453 d 
after the explosion and at a frequency of 8.4 GHz. We obtained a fairly well resolved 
image of the shell of SN 2011dh, making it one of only six recent supernovae for which 
resolved images of the ejecta are available. SN 201 ldh has a relatively clear shell mor¬ 
phology, being almost circular in outline, although there may be some asymmetry in 
brightness around the ridge. By fitting a spherical shell model directly to the visibility 
measurements we determine the angular radius of SN 2011dh’s radio emission to be 
636 ± 29 p, as. At a distance of 7.8 Mpc, this angular radius corresponds to a linear 
radius of (7.4±0.3) x 10 16 cm and an average expansion velocity since the explosion of 
I 9 OOOI 2400 km s -1 . We combine our VLBI measurements of SN 2011dh’s radius with 
values determined from the radio spectral energy distribution under the assumption 
of a synchrotron-self-absorbed spectrum, and find all the radii are consistent with a 
power-law evolution, with R ~ ^o.97±o.oi^ i m plyj ri g almost free expansion over the 
period t = 4 d to 453 d. 

Key words: supernovae: individual (SN2011dh) — radio continuum: general. 


1 INTRODUCTION 


Supernova (SN) 2011dh, discovered in the ’’Whirlpool 
Galaxy”, M51, is a recent example of a radio-loud SN. At a 
distance of 7.8Ag Mpc (Ergon et al. 2014| 1 SN 2011dh 
is also one of the nearest SNe observed in recent years. 
SN 201 ldh was discovered on 2011 May 31 by the ama¬ 


teur astronomer Amadee Riou (Griga et al. [2011) and the 
discovery was soon confirmed from pre- and post-discovery 


observations from various telescopes (Griga et al. 2011). The 


SN was coincident with the eastern spiral arm of M51. The 
explosion date is tightly constrained to be between 2011 May 
31.275 and 31.893 UT ( |Arcavi et al.|[2011| |. We will adopt 
the (rounded) midpoint of this interval, May 31.6 UT, as 
the explosion date, to, and take t to be the age of the SN 
since the explosion. 

Initially, SN 201 ldh was spectroscopically classified as 


* E-mai l: alet@hartrao .ac.za 

1 |Ergon et al.| ([2014) re-examined the various estimates of the 
distance to M51 to arrive at the value of 7.8^g'g Mpc, which is 
slightly lower than the value of 8.4 Mpc we had used in our earlier 
papers. 


Type IIP (Silverman, Filippenko & Cenko 2011), but fur¬ 


ther spectroscopy, which showed helium absorption features, 
caused a re-classification as a Type lib (Arcavi et al.|[20TT 


Marion et al. 20111. A maximum expansion velocity of 
~ 20000 km s _1 was estimated from the blue edge of the 
Ha line (Silverman, Filippenko fe Cenko|2011 Arcavi et al. 


2011 Marion et al.|20l4|. 

Radio emission was detected only a few days after the 
explosion, at centimetre wavelengths ( |Horesl i et al. 2011) 
with the National Radio Astronomy Observatory^ (NRAO) 
Karl G. Jansky Very Large Array (VLA), at millimetre 
wavelengths (Horesh, Zauderer fe Car penter|20TT I using the 
Combined Array for Research in Millimeter-wave Astron¬ 
omy, and at submillimetre wavelengths (Soderberg et al. 


20121 using the Submillimeter Array (SMA). Results from 
the initial radio and millimetre-band observations were pre¬ 
sented in Soderberg et al. (2012), while further broad-band 


measurements of the total radio flux density as well as 


2 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. 
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modelling of the lightcurve were presented in Krauss et al. 

(SSA), as is generally expected at early times (see 

Cheva- 

(2012). 

lier & Fransson 2006). Soderberg et al. (2012) and 

Krauss 

A yellow supergiant which was visible in pre-explosion 

et al. (2012) show that the radio spectrum of SN 2011dh was 


Hubble Space Telescope (HST) images, but has since disap¬ 


peared, has been identified as the progenitor (Van Dyk et al. 


2013 Ergon et al. 20141. The progenitor’s main-sequence 


mass is estimated to around 13 Mg although the various 


Anupama & Chakradhari 2013; Van Dyk et al. 2011 Maund 

et al. 2011 

Bersten et al.|2012 

). It was an extended star with 

a radius of 200 ~ 300 Ry) ( 

Van Dyk et al. 

2011 

Bersten 


et al. 2012 Ergon et al. 20141. In addition, a blue companion 
to SN 2011dh’s yellow supergiant progenitor has also likely 
been detected in deep near-UV images obtained with the 


HST (Folatelli et al. 20141, making SN 2011dh the second 


core-collapse SN, after SN 1993J (Maund et al. 2004 Fox 
et al.|2014 I, to show strong evidence of a binary companion 
for the progenitor. 

The size and expansion velocity of the shock front is 
a basic characteristic distinguishing different SNe, and it is 
therefore important to determine it observationally as di¬ 
rectly as possible. VLBI observations are the most direct 


way of making this measurement (see e.g., Marti-Vidal et al. 
2011a[|Bietenholz et al.|2010||Brunthaler et al.|2010[|Bieteir 


holz^OOli I. Unlike the optical emission, which mostly origi 

nates in the denser and more slowly moving inner ejecta, the 
radio emission generally traces the fastest ejecta. The radio 
emission is thought to originate in the region between the 
forward and reverse shocks. In the particularly well-studied 
case of SN 1993J and based on the highly resolved exem¬ 
plary shell of radio emission seen at late times, [Bartel et al.| 
12007 ) show that there is a close relationship between the 


outer boundary of the radio emission and the location of the 
forward shock at least at times > 1 yr, although |Bjornsson| 
120151 suggests that the correspondence may not be as good 


for t < 1 yr. 

In the particularly well-studied case of SN 1993J, [Bartel] 


et al. (20071 show that there is a close relationship between 


the outer boundary of the radio emission and the location 
of the forward shock. 

VLBI observations of SN 2011dh were first obtained 


at t = 14 d (2011 Jun 14, Martf-Vidal et al. 2011cI. Al 


though an accurate centre position was obtained, no useful 
constraint on the size could be obtained at this early epoch 
Further VLBI observations were obtained at epochs t = 83 
and 179 d (2011 Aug. 22 and 2011 Nov. 26) by Bietenholz 
et al. 2012 (see also [Bietenholz et al. 20TT|). At both epochs, 
we used the High Sensitivity Array (HSA), consisting of the 
NRAO VLBA, the Effelsberg telescope and the Robert C. 
Byrd Green Bank Telescope, and obtained useful measure¬ 
ments or upper limits on the source size. At t = 83 d we ob¬ 
served at 22 GHz, while at t = 179 d we observed at 8.4 GHz. 
At t = 179 d, the measured outer radius corresponds to 
an expansion velocity of 20000 ± 6500 km s -1 at 7.8 Mpc. 


Rampadarath et al. (20151 also detected SN 201 ldh in their 


1.6 GHz, wide-field VLBI observations of M51 at t = 160 d 
(2011 Nov. 7), but the source was unresolved, so again no 
useful constraint on the source radius was obtained. 

Aside from the relatively direct measurement using 
VLBI, the radius of the radio emitting region can also be de¬ 
termined from the radio spectral energy distribution (SED) 
if the spectrum is dominated by synchrotron self-absorption 


consistent with being dominated by SSA, and gave radius 
determinations up to t = 92 d. Although the calculation of 
the radius based on the SSA spectrum is fairly robust, it 
is more model-dependent than the more direct VLBI mea¬ 
surements. SN 201 ldh represents so far the best example for 
directly comparing the radii of the shock wave determined 


in these two different fashions, and Bietenholz et al. (2012) 


show that for the measurements up to t = 179 d, there is 
excellent agreement between the two methods, thereby pro¬ 
viding important confirmation for the radii derived from the 
SED by assuming SSA. 

SN 2011dh was unusual in that it remained radio-bright 
enough for VLBI observations for more than a year. We 
undertook a further epoch of VLBI observations, as well 
as observations to measure the total flux density with the 
VLA to determine the continued evolution of this SN, and 
we report on these results in this paper. 


2 OBSERVATIONS AND RESULTS 

2.1 VLA flux density observations and results 

We observed SN 2011dh with the VLA to get a total flux 
density measurements at 8.4 GHz on 2012 Aug 1 (program 
12A-286) and on 2014 Jan 31 (program 13A-370). 

The 2012 Aug 1 observations were done with the array 
in the B configuration and we used a bandwidth of 1 GHz 
centred on 8.5 GHz. They were reduced using the NRAO’s 
Astronomical Image Processing System (AIPS), with the 
flux density scale being set from observations of 3C286 us¬ 
ing the Perley-Butler 2010 coefficients. We measured an 8.5- 
GHz flux density of 0.88 ± 0.06 mjy, where the uncertainty 
includes the noise as well as an assumed 5% uncertainty in 
the flux density calibration. 

The 2014 Jan 31 observations were done with the array 
in the BnA configuration and we used a bandwidth of 1 GHz 
centred on 7.1 MHz. They were reduced using AIPS, and 
calibrated to the same flux density scale using observations 
of 3C 286. We measured a flux density of 0.66 ± 0.035 mjy, 
where again the uncertainty includes the noise as well as an 
assumed 5% uncertainty in the flux density calibration. 

In Figure [l] we show the 8.4-GHz radio lightcurve of 
SN 2011dh. We include the two new measurements described 
above as well as earlier ones taken from IBietenholz et ahl 
( |2012| ) and |Krauss et aL| (2012). The logarithmically inter¬ 
polated value at the time of our VLBI observations (2012 
Aug. 26) was 0.86 mjy. We scaled all the flux-density mea¬ 
surements taken between 7 and 9 GHz to 8.4 GHz using a 
spectral index of a? —0.7 (where S„ oc u a Kr auss et al.|2 0121. 
We are not sensitive to the exact value assumed for the spec¬ 
tral index: the earlier measurements, when the spectral in¬ 
dex may have been notably different from —0.7 were taken 
at 8.4 GHz, and therefore required no scaling, and even in 
the most extreme case of the last measurement, at 7.1 GHz, 
a difference of 0.2 in the spectral index would change the 
plotted value by less than the uncertainty. A logarithmic fit 
to the values after t — 45 d gives an average flux density de¬ 
cay rate, /?, of —0.79 ± 0.05, where S oc t~^. The lightcurve 
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Figure 1. The 8.4-GHz lightcurve of SN 201 ldh, as obtained 
from VLA observations. Flux density measurements are shown 
by solid circles, and the logarithmically interpolated value at the 
time of our VLBI observations (0.86 mjy on 2012 Aug. 26) by 
an open circle. Details for the last two measurements are given 
in the text, and the remainder are taken or interpolated from 
|Bietenholz et ah] < |2012| ) and |Krauss et al.| ( |2012| > . The age of the 
SN is calculated from the explosion date of 2011 May 31.6 UT. 


appears to flatten after t = 179 d, with a two-segment fit 
giving j3 = 1.17 ± 0.13 for t < 179 d to /3 = 0.57 ± 0.13 for 
t > 179 d. 


2.2 VLBI Observations 

We obtained VLBI imaging observations of SN 201 ldh us¬ 
ing the high-sensitivity array, which consisted of the NRAO 
Very Long Baseline Array (VLBA; 10 x 25-m diameter, dis¬ 
tributed across the United States) and the NRAO Robert 
C. Byrd (~105 m diameter) telescope and the Effelsberg 
(100 m diameter) telescope. Due to a technical failure, the 
St. Croix VLBA telescope did not observe. The observations 
occurred on 2012 Aug 26, and lasted for a total of 13 hr. At 
the midpoint of our observing session the age of the SN was 
453 d. 

To measure the tropospheric zenith delay and clock off¬ 
sets at each antenna we included three “geodetic blocks” of 
~30 minutes each at the start, middle and end of our obser¬ 
vations (see Brunthaler, Reid & Falckc 2005 Reid & Brun 


thaler 2004). In each of these geodetic blocks we observed 


9-14 bright reference sources chosen from the International 


Celestial Reference Frame (ICRF) list of sources (Fey, Gor¬ 
don & Jacobs 20091. For our geodetic blocks we recorded 


8 intermediate frequencies covering a bandwidth of 8 MHz 
spread over a ~400 MHz range. 

For the observations of SN 201 ldh we recorded a con¬ 
tiguous bandwidth of 64 MHz in each of the two senses 
of circular polarization with two-bit sampling for a total 
bit rate of 512 Mbit s . We used the same calibrator 
sources as in our previous VLBI observations of SN 201 ldh. 
We used J1332+4722 (ICRF J133245.2+472222), which is 
0.5° away from SN 2011dh, as our primary calibrator. This 
source is an ICRF source with a position known to ~70 pas 
(Fey, Gordon & Jacobs 2009). Any positions in this pa¬ 


per are calculated by taking the position of J1332+4722 to 


be Right Ascension (RA) = 13 h 32 m 45^24642, Declination 
(decl.) = 47° 22' 22'.'6670. We also obtained observations 
of the quasar, JVAS J1324+4743, about 1.5° away from 
J1332+4722, as an astrometric check source. The purpose 
of the observations of J1324+4743 was to check the quality 
of the phase-referencing and also to provide a second astro- 
metric reference source to check the positional accuracy. 

We used a cycle time of ~150 s, with ~110 s on 
SN 2011dh and ~40 s on ,11332+4722. During the observing 
run we spent three ~20 m periods observing our astrometric 
check source, J1324+4743, using a similar phase-referencing 
pattern as we used for SN 201 ldh. 

The VLBI data were correlated with the DiFX corre¬ 
lator (Deller et al.]|20111, and the analysis carried out with 
AIPS. We corrected for the dispersive ionospheric delay us¬ 
ing the AIPS task TECOR, and we solved for the zenith tro¬ 
pospheric delay on the basis of our geodetic-block observa¬ 
tions. We discarded any SN 201 ldh visibility data obtained 
when either of two of the antennas involved was observing 
at elevations below 10°. 

The initial flux density calibration was done through 
measurements of the system temperature at each telescope, 
and then improved through self-calibration of the primary 
reference source J1332+4722. This source is slightly re¬ 
solved, as can be seen on the images in the VLBA calibrator 
list data-basB and also from our previous observations at 
22 and 8.4 GHz (Bietenholz et al. 20121. We see a simi¬ 
lar structure in the image from our current observations at 
8.4 GHz, where a weak extension or second component is 
visible ~2 mas to the west-southwest of the peak. 

Our final amplitude and phase calibration was de¬ 
rived using a CLEAN model of this source, with the peak- 
brightness point in the image being placed at the nominal 
coordinates given abo\ r B Finally this calibration was inter¬ 
polated to the intervening scans of SN 201 ldh. 


2.3 VLBI Results 


We detected SN 201 ldh with a sufficient signal-to-noise ra¬ 
tio to obtain a high-quality VLBI image with a resolution 
of 0.79 x 0.52 mas (FWHM), a dynamic range of 23, and an 
image background rms of 12 /xjy beanU 1 . We show the im¬ 
age in Figure[2] The total CLEAN flux density was 760 /rJy, 
which was ~90% of the total flux density measured using 


the VLA (see section 2.1l. The structure of SN 2011dh is 


relatively circular at the lower contours, but two hot-spots, 
bilaterally located approximately east and west, are evident 
in the image. 


As in our previous VLBI observations (Bietenholz et al. 


2011 2012), we again turn to fitting a spherical-shell model, 


with an outer radius of 1.25X the inner radius, directly to 
the visibility measurements to accurately measure the size 
and centre position of the source. The free parameters in the 
fit are the centre position, the radius, and the flux density. 

We take the fitted centre position of the model as our 


3 http://www.vlba.nrao.edu/astro/calib 

4 The deviation of the source geometry of J1332+4722 from a 
point source is small enough so that the effect of using a point 
model in the solutions for delay and delay rate made using FRING 
is negligible. 
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Figure 2. The 8.4-GHz VLBI CLEAN image of SN 2011dh on 
2012 August 26, at t = 453 d. The greyscale and contours both 
show the brightness. The greyscale is labelled in ^zjy beam -1 , and 
the contours are drawn at —10, 10, 20, 30, 50 (emphasized), 70, 
80 and 90% of the peak brightness, which was 270 ^Jy beam -1 . 
The convolving beam, indicated at lower left, was 0.79 X 0.52 mas 
(FWHM) at p.a. —15°. The image was made using Briggs’ robust 
weighting ( jBriggs, Schwab Sz Sramek| 1 999| , with the weight¬ 
ing set to approximately midway between uniform and natu¬ 
ral weighting (AIPS robustness parameter set to —1). The to¬ 
tal CLEAN flux density was 760 [i Jy, and the background rms 
level was 12 fijy beam -1 . North is up and east is to the left, and 
the coordinates are given as offsets from the centre position of 
the fit model, which is at RA = 13 h 30 m 5? 1055478 and decl. = 
47° 10' 10'.'92261 (see Text.[Z3]i. 



best estimate of the centre position of SN 201 ldh. For the 
observations of SN 2011dh, at t = 453 d, the position of 
the model centre was RA = 13 h 30 m 5U055478 and decl. = 
47° 10 7 10 , /92261. The statistical uncertainty on this posi¬ 
tion is 30 /ras in both RA and decl. 

The fitted outer radius was 636 ± 29 /xas. The statis¬ 
tical uncertainty was ±23 /ras, and a Monte-Carlo simula¬ 
tion with 16 trials where we randomly varied the amplitude 
gains by a = 20% gave an additional radius uncertainty due 
to possible amplitude calibration errors of 13 /ras. We add 
these contributions in quadrature, and our final standard 
error on the fitted outer angular radius is therefore 29 fias. 
At 7.8lj g Mpc, this angular radius corresponds to a linear 
radius of (7.4lQ g) x 10 16 cm, and an average expansion ve¬ 
locity since the explosion of I 9 OOOI 2400 km s -1 , where our 
uncertainties on the linear radius and expansion velocity in¬ 
cludes the contribution of the uncertainty of the distance. 


3 DISCUSSION 

We obtained phased-referenced VLBI observations of 
SN 2011dh, with the primary goal of obtaining a clearly 
resolved image of the expanding shell of ejecta, at t = 453 d 


after the explosion. We obtained a fairly well resolved im¬ 
age of SN 2011dh, shown in Fig [2] which shows a shell like 
structure that is fairly circular. Our image shows a bilateral 
enhancement of brightness, with hot-spots located approxi¬ 
mately east and west along the ridge. 

In the case of the similar Type lib SN, SN 1993J, 
there are clear time-dependent modulations of the bright¬ 
ness along the ridge-line which were larger than the noise or 


any systematic effect during the first two years ( Bietenholz, 
|Bartel &; Rupen|2003||Bietenholz|2008[ ). Do the hot-spots in 
the image of SN 2011dh indicate a real brightness enhance¬ 
ment along the ridge? 

Even in the case of a completely circular structure, such 
a bilateral enhancement is in fact expected in the case of an 
elliptical CLEAN beam, with the hot-spots occurring at a 
p.a. at right angles to the elongation of the beam. We made 
simulated visibility data sets of a purely circular shell model 
with a noise level similar to that in the real data. When 
these were CLEANed in a fashion similar to the real data, 
the resulting images also tended to display apparent hot 
spots, similar in brightness to those visible in Figure [2] We 
therefore conclude that although the bilateral enhancement 
seen in Figure [2] could be real, a perfectly uniformly circular 
structure is also compatible with the data. An enhancement 
on the west side is suggested, but not demanded by the data. 

Optical polarization measurements of SN 201 ldh sug¬ 
gest significant departures from sphericity for times t < 30 d 
(Mauerhan et al.j20151, which is in contrast to the relatively 
circular outline seen in our VLBI image. A similar contrast 
between a circularly symmetric VLBI image and optical po¬ 
larization suggesting significant departures from symmetry 
exists in SN 1993J. The reason for the contrast is likely that 
the radio emission arises from the outer shock, which is rel¬ 
atively spherical, while the optical polarization arises from 
the inner ejecta, which are much less so. 

We combined the present VLBI measurement of the 
centre position and outer radius with those from earlier ob¬ 
servations from Bietenholz et al. (20121, Krauss et al. (20121 
and Martf-Vidal et al. (2011cI. 

Since slightly different correlator positions for the 
phase-reference source J1322±4722 were used in the differ¬ 
ent VLBI experiments, we corrected the phase-referenced 
SN 2011dh centre position as well as that of the check source, 
J1324±4743, so that they are all relative to the position of 
J1322±4722 given in §2.2| above. 

Since all three of our VLBI runs included observations 
of the check source, we can compare the phase-referenced 
position of J1324±4743 among the three epochs. We found 
that the scatter in this position over the three epochs was 
130 fias in R.A. and 90 /ras in decl. This scatter is larger than 
the expected astrometric errors, suggesting that there is 
some small apparent motion in one (or both) of the sources. 
Since we do not know whether our phase-reference source, 
J1324±4743, or the check source J1322±4722 has an appar¬ 
ently variable position, we take this scatter to be the uncer¬ 
tainty in the position of the phase-reference source. Based 
on this, we then take a conservative value of 130 /ias for the 
uncertainty also in the centre position of SN 2011dh. 

Turning now to the centre positions of SN 2011dh, and 
including the value at t = 14 d from |Martf-Vidal et al.| 
(2011c I, we performed a least-squares fit of the position as a 
function of time. We find a proper motion of 55±66 fj ,as yr _1 
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Figure 3. The shock front radii of SN 201 ldh as estimated by two 
independent methods. Using red circles, we plot the values derived 
from fitting spherical shell models to the VLBI visibility-data 
from this paper and |Bietenholz et al.| ( (2012[ ). The plotted lcr er¬ 
ror bars include statistical and systematic contributions (see Sec¬ 
tion [53] for details) but do not include the contribution from the 
distance, since any error in the distance would affect all the points 
in the same way. It should be noted that the error bars on the 
last point are smaller than the point symbol and is therefore not 
visible in the plot. Using triangles, we plot the values calculated 
from the radio spectral energy distribution under the assumption 
that it is dominated by SSA (synchrotron self-absorption), with 
green triangles pointing up for the values taken from [Soderberg] 
|et al.| J 20 T 2 I “S+2012”) and |Krauss et al.| ( [20121 “K+2012”) and 
blue triangles pointing down for those from |Horesh et ah ( 2013 J. 
The uncertainties include statistical and systematic components 
and all radius measurements have been rescaled to our distance 
of 7.8 Mpc. The blue line represents our power-law fit to all the 
values, with R = 5.30 X 10 15 (t/30d)° 97 cm (see text). 


(at p.a. —15°), which corresponds to 2000 ± 2400 km s -1 , 
with a 3 <t upper limit of 9200 km s -1 . 

We note that any proper motion due M51’s galactic 
rotation is expected to be only a few hundred km s _1 , well 
below our uncertainty, so we can take this projected speed 
to be that of SN 2011dh’s centre with respect to its local 
frame of rest within M51. 

For the outer radius of the radio emitting region, which 
as we have argued, should be closely tied to the outer shock 
radius, we use the value determined in the present work 
along with the two VLBI values from |Bietenholz et al~| 
(20121, and those determined from the radio spectrum by 
assuming SSA from Krauss et ah] (20121, as as well as from 


Horesh et al. (20131 (t ~ 4toll d). We scale all the radii to 


a consistent distance of 7.8 Mpc. We plot them in Figure [3] 
The expansion of a SN is usually parametrized as a 
power-law with R = At" 1 , where R is the radius, t is the 
time, m is the power-law index, called the “deceleration pa¬ 
rameter” and A is the radius at t = 1. A weighted least- 
squares fit to all the values of R gives 

R = (5.3 ± 0.1) (t/30 d) (0 ' 97±0 01) x 10 15 cm 


The uncertainty on m is the statistical one only, calcu¬ 
lated assuming that the points are all independent. All the 
radius measurements are consistent with a power-law evo¬ 
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lution. The fitted expansion index m = 0.97 ± 0.01 implies 
almost free expansion. 

As already noted in Bietenholz et al. (21) 12), the radii 
measured with VLBI are consistent with those determined 
somewhat less directly from the SED under the assump¬ 
tion of SSA. We further note that the SSA radii derived 


by Horesh et al. (2013) are quite consistent with those of 


Soderberg et al. (2012) and Krauss et al. (2012). The ear¬ 


lier VLBI radius measurements were rather more uncertain 
than those obtained from the SED. Our new measurement, 
by contrast, is rather more precise, and is still consistent 
with the powerlaw implied by the SED-based radius mea¬ 
surements. Note that at t > 100 d, the turnover frequency 
becomes too low to be directly measured, thus it is not pos¬ 
sible to obtain radii from SED fits at the times of our more 
accurate VLBI measurements. Nonetheless, while it is possi¬ 
ble that the SED-derived radii are systematically too small 
and that the expansion decelerated after t ~ 100 d, we con¬ 
sider it an unlikely coincidence that the measured VLBI ra¬ 
dius at t = 453 d would lie so close to the extrapolation of 
the powerlaw seen at early times in that case. 

The radii determined from the SED depend on D (18 / 19 ^ 
whereas the (linear) radii determined from the VLBI mea¬ 
surements depend on D, so only very large errors in D would 
produce a substantial scaling between the two kinds of ra¬ 
dius determinations. We note, however, that by comparing 
the angular expansion velocity as determined from VLBI 
and less directly from the SED fits, with linear expansion 
velocity measured from optical spectroscopy, a direct dis¬ 
tance can be determined via the “expanding shock front” 


method. See Bartel et al. (2007) for a description of this 


method, and its application in the case of SN 1993J in M81. 
We defer a detailed comparison along these lines to a future 
paper, but note here that our fit implies an expansion ve¬ 
locity at t = 30 d of 2 OOOOI 2400 km s_1 j where again our 
uncertainty on the expansion velocity includes the contri¬ 
bution from the distance uncertainty (+15%, —12%). This 
velocity compares very well to maximum velocity seen at the 
blue edge of Ho absorption seen in the optical at the same 
age (e.g. Marion et al. |201 4). 


3.1 Self-similar approximation and estimate of 
CSM and ejecta density profiles 

The amount of radio emission as well as the deceleration are 
ultimately dependent on the energy dissipated as the shock 
ploughs out through the CSM, and depend therefore on the 
density profiles in the ejecta and the CSM. In particular, if 
the density profiles of the ejecta and the CSM are power- 
laws in radius, then a self-similar solution exists, and both 
the total flux-density and the shock radius follow powerlaws 
in time such that the flux density decays as S oc while 
r oc t m (Chevalier 1982 Fransson, Lundqvist & Chevalier 


1996). 


Our fit to the expansion of SN 201 ldh gives m = 0.97± 
0.01, and our fit to the optically thin part of the 8.4-GHz 
lightcurve gives /? = —0.79 ± 0.05. If we assume power-law 
density profiles, with the ejecta density being oc r" and the 
CSM density oc r s , and taking the radio spectral index to 
be a = —0.7 ± 0.2, then those values of m and /3 imply 
n = 31 ± 11 and s = 2.0 ± 0.1. So the implied CSM density 
profile overall is consistent with a wind density profile with 
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p oc r~ . The ejecta density profile is less well determined, 
but the small deceleration and flat lightcurve imply a very 
steep profile in the ejecta, consistent with what was found 


by Maeda (20120 based on X-ray observations. 


Although our radius measurements are compatible with 
an unbroken powerlaw, the measurements suggest a flat¬ 
tening after t — 179 d instead of a single powerlaw, with 
p = 1.17T0.13 until then and 0.57±0.13 thereafter. No ob¬ 
vious corresponding change is seen in the expansion curve. 
A flattening of the flux density decay implies either a steeper 
density profile in the ejecta, or a flatter one in the CSM. The 
former would be accompanied by a decrease in deceleration 
while the latter would be accompanied by an increase. 

Since we see no obvious change in the deceleration at t = 
179 d, our measurements do not reliably distinguish between 
these two cases. Although the average deceleration over all 
the measurements is well determined atm = 0.97±0.01, the 
value for t > 179 d is much less so, and a change of up to 
0.1 in m at t > 179 d is compatible with our measurements. 
Since m > 1 can only happen with some form of continued 
energy input, we think it more likely that m has decreased 
at the latest times. 

We note that a density structure more complicated than 
a simple powerlaw is not unreasonable for either the ejecta 
nor the CSM, and that in fact departures from powerlaw 
density structure have been suggested for both the ejecta 
and the CSM of the well-studied Type lib SN 1993J. We 
will elaborate on the comparison between SN 201 ldh and 


that such local increases in the ejecta density can cause 
temporary flattening of the lightcurve decay. |Bartel et al.| 
12002) looked at the radio lightcurve and expansion rate of 


SN 1993J till t = 3164 d and suggest changes in the CSM 
density profile. 

Turning to SN 2011dh, overall the fast rise of the radio 
lightcurve to its maximum, the small deceleration and the 
low peak luminosity imply that the CSM in this case was 
less dense than that of SN 1993J. Other Type lib SNe which 
have shown similarly quick rise times and relatively low ra¬ 
dio luminosities are SN 2008bo ( Stockdale et al.||2008 l, and 
SN 2008ax (Romin g et al.|2009 l, and indeed the similarities 
of SN 2011dh’s optical spectra to those of SN 2008bo have 
already been pointed out (Maund et al.|2011 Marion et al. 
20141 


The small deceleration suggests that the swept-up mass 
must still be low compared to that of the ejecta. Kraus set al.| 
(20121, Horesh et al. (20131, and Maeda et al. (20151 have all 


estimated mass-loss rates for SN 2011dh’s progenitor, with 
the largest of those three estimates being that of Maeda et al, 
which was ~ 3 x 10 -5 Mg yr -1 . Using that mass-loss rate, 
and a wind speed for the yellow supergiant of 20 km s -1 , 
we calculate that the swept-up mass at t = 453 d was 
1.3 x 10 -3 Mq. If we used the mass-loss rates of 


Krauss et al. 


SN 1993J in the following section (13.21. 


( |2012[ ) or |Horesh et al.| ( |2013| we would derive a somewhat 
lower swept-up mass. Since the ejected mass is still low with 
the mass of the ejecta (a few Mg) at t = 453 d, our finding 
that the SN is still in almost free expansion is consistent 
with the expectations. 


3.2 Comparing SN 2011dh to SN 1993J. 

We can compare SN 2011dh to SN 1993J: both were of 
Type lib, and had extended progenitors with binary com¬ 
panions, and had lost a considerable fraction of their hydro¬ 
gen envelopes to the companion before the SN explosion. For 
SN 1993J, the radio lightcurves and expansion were partic¬ 


ularly well determined (see, e.g. Bartel et al. 2002 Martf- 


Vidal et al. 2011b). SN 1993J was about four times as ra¬ 


dio luminous as SN 201 ldh, as the former reached a peak 
8.4-GHz spectral luminosity of ~ 2 x 10 * 2 3 4 5 ' erg s -1 Hz -1 , 
compared to only 0.5 x 10 2 ' erg s _1 Hz -1 for the latter. 
The evolution of SN 1993J’s radio lightcurves was much 
slower than that of SN 2011dh, since the former reached 
an 8.4 GHz peak only after t ~ 150 d, compared to the 
t = 30 d for SN 201 ldh. However, we note that a circum- 
stellar density profile flatter with s < 2, i.e. flatter than that 


of a steady wind, was also suggested for SN 1993J (e.g. van 


|Dyk et al.|1994||Fransson, Lundqvist fc Chevalier|1996| ), al¬ 

though models with a constant s = 2 can also reproduce the 


measurements (e.g. Fransson & Bjornsson 1998 Martf-Vidal 
|et al.|2011b| ). 

In the case of SN 1993J, departures from simple pow¬ 
erlaw profiles have been suggested for both the ejecta and 
the CSM. Iwamoto et al. (19971 modelled the optical emis¬ 


sion, while Suzuki & Nomoto (1995) modelled the early X 


ray emission of SN 1993J, and both authors found that to 
reproduce the observed emission required an ejecta struc¬ 
ture considerably more complicated than a simple power- 
law, including an increase in density going inward past the 
H/He interface. Modelling of the radio lightcurves and ex¬ 
pansion by Mioduszewski, Dwarkadas & Ball (2001) showed 


4 CONCLUSIONS 

1. We obtained a new, resolved, VLBI image of SN 2011dh, 
at t = 453 d after the explosion. We found that SN 2011dh 
shows a relatively circular projected shell structure, possi¬ 
bly with some enhancement on the west side. With our new 
VLBI image, SN 201 ldh takes its place among the only six 
recent SNe for which resolved images of the ejecta are avail¬ 
able. With our new VLBI image, SN 201 ldh takes its place 
among the only six recent SNe for which resolved images 
of the ejecta are available — the others being SN 1979C 


(Bartel & Bietenholz 20081, SN 1987A (Ng et al. 20111, 


SN 1986J (Bietenholz, Bartel & Rupen 20101, SN 1993J 
( |Bietenholz, Bartel fe Rupen||2003| |Marcaide et al.|[l997| ), 

SN 2008iz ( Brunthaler et al.|2010 i. 

2. We also obtained measurements of the total flux density 
with the VLA, and the radio lightcurve shows an approx¬ 
imately power-law decay with S oc t -0,79 ), although the 
decay flattens somewhat after t = 179 d. 

3. Our astrometric VLBI measurements give an upper limit 
of 9200 km s _1 on the average projected speed of the cen- 
trepoint of SN 2011dh since the explosion. 

4. We measured a radius of SN 2011dh at t m 453 d of 
(7.4±S;J) X 10 16 cm using our VLBI data. Comparing to 
earlier values derived both from VLBI observations and from 
fitting of the SED, we show that SN 2011dh’s expansion is 
well fit with a power-law with R oc ^°- 97±01 °. The SN is 
therefore still in almost free expansion. 

5. Our fit to the expansion implies a velocity for the forward 
shock at t = 30 d of 20 0 00))24oo km s _1 , which is comparable 
to the highest velocities seen in optical spectral lines. 
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6. The expansion curve and the radio lightcurve together are 
consistent with a wind density profile (p oc r -2 ) in the cir- 
cumstellar medium, but suggest a very steep density profile 
in the ejecta. 
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